Ab initio calculations of the dynamical response of copper 
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The role of localized d-bands in the dynamical response of Cu is investigated, on the basis of ab initio 
pseudopotential calculations. The density-response function is evaluated in both the random-phase 
approximation (RPA) and a time-dependent local-density functional approximation (TDLDA). Our 
results indicate that in addition to providing a polarizable background which lowers the free-electron 
plasma frequency, d-electrons are responsible, at higher energies and small momenta, for a double- 
peak structure in the dynamical structure factor. These results are in agreement with the exper- 
imentally determined optical response of copper. We also analyze the dependence of dynamical 
scattering cross sections on the momentum transfer. 

PACS numbers: 71.45. Gm, 72.30.+q, 78.20.-e, 78.70.Ck 



Noble-metal systems have been the focus of much ex- 
perimental and theoretical work in order to get a bet- 
ter understanding of how electronic properties of delocal- 
ized, free-electron-like, electrons are altered by the pres- 
ence of localized <i-electrons. Silver is one of the best 
understood systems, where the free-electron plasma fre- 
quency is strongly renormalized (red-shifted) by the pres- 
ence of a polarizable background of d-electronstl. Un- 
like silver, copper presents no decoupling between sp 
and d orbitals, and a combined description of these one- 
electron states is needed to address both structural and 
electronic properties of this material. Though in the 
case of other materials, as semiconductors, electron-hole 
interactions (excitonic renormalization) strongly modify 
the single-particle optical absorption profileo, metals of- 
fer a valuable playground for investigations of dynamical 
exchange-correlation effects of interacting many-electron 
systems within a quasiparticle picturea. Indeed, ab ini- 
tio calculations of the dynamical response of a variety 
of simple metals, as obtained within the random-phasc- 
approximation (RPA), successfully account for the exper- 
imentally determined pLasmon dispersion relations^ and 
scattering cross sectionaj. Within the same many-body 
framework, ab initio calculations of the electronic stop- 
ping power of real solids have also been reportedH. 

Since_the pioneering investigations by Ehrenreich and 
PhilippQ on the optical absorption and reflectivity of Ag 
an Cu, a variety of measurements of the optical properties 
of copper has been reportedQ. Nevertheless, there have 
been, to the best of our knowledge, no first-principles 
calculations of the dielectric response function of Cu that 
include the full effects of the crystal lattice. Furthermore, 
the dynamical density-response function is well-known to 
be a key quantity in discussing one-electron properties 
in real metals, and it is also of basic importance in the 
description of low-dimensional copper systems studieeLin 
optical and time-resolved femtosecond experimentsoEJ. 

In this Rapid Communication we report a first- 



principles evaluation of the dynamical density-response 
function of Cu, as computed in the RPA and a time- 
dependent extension of local-density functional theory 
(TDLDA), after an expansion of 4s 1 and 3c? 10 one- 
electron Bloch states in-a plane- wave basis with a kinetic- 
energy cutoff of 75 RytiJ. Though all-electron mixed basis 
schemes, such as the full-potential linearized augmented 
plane wave (LAPW) methodEj, are expected to be well 
suited for the description of the response of localized d- 
electrons, plane-wave pseudopotential approaches offer a 
simple and well-defined scenario to describe gxound-state 
properties and dynamical response functionsEj. This ap- 
proach has already been successfully incorporated in the 
description of inelastic lifetimes of excited electrons in 
coppeiO, and could also be applied to the study of other 
noble and transition metals. 

The key quantity in our calculations is the dynam- 
ical density-response function %(r, r';w). For periodic 
crystals we Fourier transform this function into a matrix 
Xg,G'{<\, which, in the framework of time-dependent 
density-functional theory (DFT)E£rE3, satisfies the matrix 
equation 



ag,g' (q, uj) = xg,g> (q, w ) + mi ag,g" (q, w) 

G" G'" 

X [v G "(q)<$G'',G'" + #G»,G"'(q. W )] XG"',G'(q,w). (1) 



Here, the wave vector q is in the first Brillouin zone (BZ), 
G and G' are reciprocal lattice vectors, «G(q) = 47r/|q+ 
G| 2 are the Fourier coefficients of the bare Coulomb po- 
tential, the kernel A' G c G ,(q j w) accounts for short-range 
exchange-correlation effectsEJ, and x GG ,(q, lj) are the 
Fourier coefficients of the density-response function of 
noninteracting Kohn-Sham electrons: 
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where the second sum runs over the band structure for 
each wave vector k in the first BZ- f^ n are Fermi fac- 
tors, r\ is a positive infinitesimall3, and represents 
the normalization volume. The one-particle Bloch states 
<^k,n(r> to) and energies i?k,n are the self-consistent eigen- 
functions and eigenvalues of the Kohn-Sham equations of 
DFTlia, w hich we solv£_within the so-called local-density 
approximation (LDA)E3 for exchange and correlation ef- 
fects. The electron-ion interaction is described-by a non- 
local, norm-conserving, ionic pseudopotentiaEj in a non- 
separable form, so as to have a better description of the 
conduction bands entering Eq. (||). 

The calculations presented below have been found to 
be well converged for all frequencies and wave vectors 
under study, and they have been performed by includ- 
ing conduction bands up to a maximum energy of 80 eV 
above the Fermi level. BZ integrations were performed by 
sampling on a 10 x 10 x 10 Monkhorst-Pack meshta In the 
RPA the kernel i"T G c G ,(q, to) is taken to be zero. In the 
TDLDA the zero-frequency kernel, approximated within 
the LDA by a contact delta iupction, is adiabatically 
extended to finite frequenciesliSEi In both RPA and 
TDLDA, crystalline local field effects appearing through 
the dependence of the diagonal elements of the interact- 
ing response matrix XG,c(q,w) on the off-diagonal el- 
ements of the polarizability Yq g , (q, ui) have been fully 
included in our calculations^!. 

The properties of the long-wavelength limit (q — > 0) 
of the dynamical density-response function are accessi- 
ble by measurements of the optical absorption, through 
the imaginary part of the dielectric response function 
£G=o,G'=o(q = 0, to). On the other hand, the scatter- 
ing cross-section for inelastic scattering of either X-rays 
or fast electrons with finite momentum transfer q + G is, 
within the first Born approximation, proportional to the 
dynamical structure factor 



S(q + G,tu) 
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where 



e G 1 ,G'(^ UJ ) = S G,G> +«G'(q)XG,G'(q,w). 
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Fig. 1 exhibits, by solid lines, our results for both real 
and imaginary parts of the eG,G(q, to) dielectric function 
of copper, for a small momentum transfer of |q + G| = 
0.18 clq 1 (ao is the Bohr radius), together with the optical 
data (q = 0) of Ref. || (dashed lines). In this low-q limit, 
both RPA and TDLDA dynamical density-response func- 
tions coincide, and the dielectric function is obtained 
from the dynamical density-response function of nonin- 
teracting Kohn-Sham electronsa. Corresponding values 

of the so-called energy-loss function Im — e G 1 G (q, to) 

are presented in Fig. 2, and a comparison between the 
imaginary parts of interacting xg.g(cIiw) and noninter- 
acting Xg g (q> u ) density-response functions is displayed 
in the inset of this figure, showing that as the Coulomb 



interaction is turned on the oscillator strength is redis- 
tributed. Our results, as obtained for a small but fi- 
nite momentum transfer, are in excellent agreement with 
the experimentally determined dielectric function, both 



£ G.G 



(q,w) 



showing a double peak structure in Im 

In order to investigate the role of localized d-bands 
in the dynamical response of copper, we have also used 
an ab initio pseudopotential with the 3d shell assigned 
to the core. The result of this calculation, displayed 
in Fig. 2 by a dotted line, shows that a combined de- 
scription of both localized 3d 10 and delocalized 4s 1 elec- 
trons is needed to address the actual electronic response 
of copper. On the one hand, the role played in the 
long-wavelength limit by the Cu d-bands is to provide 
a polarizable background which, lowers the free-electron 
plasma frequency by ~ 2.5eVc3. We note from Fig. 1 
that near 8.5 eV the real part of the dielectric function 
(Ree) is zero; however, the imaginary part (Ime) is not 
small, due to the existence of interband transitions at 
these energies which completely damp the free-electron 
plasmon. On the other hand, d-bands are also responsi- 
ble, at higher energies, for a double-peak structure in the 
energy-loss function, which stems from a combination of 
band-structure effects and the building up of collective 
modes of d-electrons. Since these peaks occur at energies 
(~ 20 eV and ~ 30 eV) where Re e is nearly zero (see Fig. 
1), they are in the nature of collective excitations, the 
small but finite value of Im e at these energies accounting 
for the width of the peaks. 

A better insight onto the origin of the double-hump in 
the energy-loss function is achieved from Fig. 3, where 
the density of states (DOS) and the joint-density of states 
( J-DOS) of Cu are plotted. The high-energy peak present 
in the J-DOS spectrum at about 25 eV, which appears as 
a result of transitions between d-bands at ~ 2 eV be- 
low the Fermi level and unoccupied states with energies 
of ~ 23 eV above the Fermi level, is responsible for the 
peak of electron-hole excitations in Ime and Imx° at 
lo = 25 eV (see Fig. 1 and the inset of Fig. 2). Hence, 
there is a combination, at high energies, of d-like collec- 
tive excitations and interband electron- hole transitions, 
which results in a prominent double-peak in the loss spec- 
trum. 

Now we focus on the dependence of the energy-loss 
function on the momentum transfer q+G. As long as 
the 3d shell of Cu is assigned to the core, we find a well- 
defined free-electron plasmon for wave vectors up to the 
critical momentum transfer where the plasmon excitation 
enters the continuum of intraband particle-hole excita- 
tions. This free-electron plasmon, which shows a charac- 
teristic positive dispersion with wave vector, is found to 
be completely damped when a realistic description of 3d 
orbitals is included in the calculations. At higher energies 
and small momenta, d-like collective excitations originate 
a double-peak structure which presents no dispersion, as 
shown in Fig. 4. In this figure the RPA dynamical struc- 
ture factor for G = is displayed, as obtained for various 
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values of q along the (100) direction. For larger values of 
the momentum transfer q+G, single-particle excitations 
take over the collective ones up to the point that above a 
given cutoff the spsGtea is completely dominated by the 
kinetic-energy termcHI9. 

In Fig. 5 we show the computed dynamical struc- 
ture factor for |q + G| = l.Etla^ 1 along the (111) di- 
rection, in both RPA (solid line) and TDLDA (dashed 
line), together with the result of replacing the interact- 
ing xg,g(si,u) matrix by its noninteracting counterpart 
XQ G (q, uj) (dotted line). The noninteracting dynamical 
structure factor (dotted line) now reproduces the main 
features of full RPA and TDLDA calculations. The dou- 
ble peak of Fig. 2, which is in the nature of plasmons, 
is now replaced by a less pronounced double-hump orig- 
inated from single electron-hole excitations. A similar 
double-peak has been found in the loss spectra of sim- 
ple metalsEl, which has been understood on the basis_af 
the existence of a gap region for interband transitionscj. 
We also note that the effect of short-range exchange- 
correlation effects, which are absent in the RPA, is to 
reduce the effective electron-electron interaction, thus the 
dynamical structure factor being closer in TDLDA than 
in RPA from the result obtained for noninteracting Kohn- 
Sham electrons. The RPA dynamical structure factor of 
Cu is enhanced by up to a 40% by the inclusion, within 
the TDLDA, of many-body local field corrections. 

In summary, we have presented ab initio pseudopoten- 
tial calculations of the dynamical density-response func- 
tion of Cu, by including d-electrons as part of the valence 
complex. In the long- wavelength limit (q — > 0), d-bands 
provide a polarizable background that lowers the free- 
electron plasma frequency, d-electrons are also responsi- 
ble for a full damping of this s-like collective excitation 
and for the appearance of a d-like double-peak structure 
in the energy-loss function, in agreement with the experi- 
mentally determined optical response of copper. We have 
analyzed the dependence of the dynamical structure fac- 
tor on the momentum transfer, and we have found that, 
for values of the momentum transfer over the cutoff wave 
vector for which collective excitations enter the contin- 
uum of intraband electron-hole pairs, a less-pronounced 
double-hump is originated by the existence of interband 
electron-hole excitations. Experimental measurements of 
scattering cross sections in Cu would be desirable for the 
investigation of many-body effects, which we have ap- 
proximated within RPA and TDLDA. 
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FIG. 1. Real and imaginary parts of the ea,G,(ci,u}) dielec- 
tric function of Cu, for q = 0.2(1,0,0) x (2ir)/a (a = 3.61 A 
is the experimental lattice constant) and G = 0. Solid and 
dashed lines represent our calculations and the optical data 
of Ref. [], respectively. 

FIG. 2. As in Fig. 1, for the energy-loss function of Cu. 
The dotted line represents the result of assigning the 3d shell 
to the core. In the inset the imaginary parts of interacting 
(solid line) and noninteracting (dashed line) density-response 
functions are compared, for the same values of q and G. 

FIG. 3. Calculated Density of States (DOS) and 
Joint-Density of States ( J-DOS) of Cu. The 25 eV peak in the 
J-DOS mainly corresponds to transitions from the d band to 
unoccupied states with energies of ~ 23 eV above the Fermi 
level. The J-DOS is a first order approximation to the optical 
spectra, when matrix-element renormalization and Coulomb 
interactions are neglected. The DOS and the J-DOS were 
computed using an interpolataion scheme based on the linear 
tetrahedron method with a 20 x 20 Mokhrost-Pack mesh. 

FIG. 4. RPA energy-loss function of Cu along the (100) 
direction, for various values of the momentum transfer |q+G|: 
0.2, 0.4 and 0.8, in units of 2-w/a (a = 3.6lA). 

FIG. 5. RPA (solid line) and TDLDA (dashed line) dy- 
namical structure factors of copper, for |q + G| = l.Qla^ 1 
along the (111) direction (q = 0.2(1,1,1) x (2ir/a) and 
G = (1, 1, 1) x (2n/a)). The dotted line represents the re- 
sult of replacing the interacting density-response matrix by 
its noninteracting counterpart. 
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